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TECHNICAL NOTE 2375 


ON THE USE OF COUPLED MODAL FUNCTIONS .IN FLUTTER ANALYSIS 
By Donald S. Woolston and Harry L. Runyan 


SUMMARY 


An investigation of the flutter characteristics of a uniform, 
unswept, cantilever wing .of high aspect ratio and under conditions of 
high mass coupling has been made by means of an analysis of the Rayleigh 
type based on coupled modal functions. Results are compared with experi- 
ment and also with the calculated results of NACA TN 1902 in which 
uncoupled modes were used. For the configuration studied, the use of 
coupled modes yielded, in general, no better agreement with experiment 
than did the use of uncoupled modes and, in some cases, the uncoupled- 
mode approach was better. 


INTRODUCTION 


In the study of flutter, many simplifications must be made in 
order to obtain practical solutions. One simplification common to most 
generally used methods of analysis is the use of a finite series of 
modal functions to represent wing motion during flutter. Of the many 
functions which could be used, either the coupled or the uncoupled modes 
of vibration of the system are usually chosen. A question, which natu- 
rally arises and one which has been a matter of practical interest to 
flutter analysts for some time, is that of the better choice of these 
two approaches. The coupled-mode approach has been found to be very 
time consuming compared with the uncoupled-mode approach. Various 
investigators (see, for example, reference l), however, have expressed 
the thought that the use of coupled modes might give a higher degree of 
accuracy that would compensate for the greater amount of labor involved. 

The present paper, which deals with the use of coupled modes, and 
references 2, 3 , and 4 investigate this question. Reference 2 presents 
the experimental results of an extensive testing program made with an 
unswept, uniform, cantilever wing of fairly high aspect ratio. Flutter 
tests were made with a single concentrated weight mounted at various 
spanwise and chordwise positions on the wing. 

In reference 3 a differential-equation analysis was applied to 
some of those cases of reference 2 where large mass coupling was 
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involved. Good agreement between theory and experiment was obtained 
for all cases studied. The results indicated that the structural part 
of the problem was adequately taken into account by the differential- 
equation approach and that the theoretical two-dimensional aerodynamic 
coefficients were also adequate for the conditions investigated. 

In reference 4 an analysis of the Rayleigh type, based on 
uncouple d-mo dal functions, was applied to a number of the experimental , 
cases of reference 2. An effort was made to appraise the accuracy of 
this method of analysis and to determine the number of uncoupled modes 
needed to give a satisfactory result. No guide as to the number of 
modes which should be considered could be given and, in some cases, 
more than a practical number seemed necessary. 

The present investigation is an extension of the work of these 
references and deals with the application of a Rayleigh type of analysis 
in which coupled modes are used. Analyses have been made for several 
spanwise stations of the case in reference 2 which is designated by 
weight 7 and leading-edge position a, as well as for the wing without 
a weight. Weight "Ja was selected for analysis since, for this partic- 
ular weight, the least satisfactory agreement with experiment was 
obtained by the uncoupled-mode approach (reference 4) . The results of 
the present investigation in which coupled modes are used are compared 
with the results of the uncoupled-mode approach of reference 4. Since 
the results of reference 3 indicated that two-dimensional aerodynamic 
coefficients were adequate for these cases, this direct comparison of 
modal approximations seems possible. 

A general outline of the procedure involved in conducting a 
coupled-mode analysis is given and includes the form of the flutter 
determinant . Application of the method to the specific cases is then 
discussed. 


SYMBOLS 


a nondimensional distance of elastic axis from midchord 

measured in half-chords, positive for positions of elastic 
axis behind midchord 

A^ flutter- determinant element associated with kinetic and 

potential energies of mechanical system 

b wing half-chord 

Cj^ flutter- determinant element associated with energies of air 

stream 



flutter frequency, cycles per second 

structural damping coefficient considered as variable in 
solution of flutter determinant 

structural damping coefficient in jth coupled mode 

bending component of ith coupled mode of vibration 

mass moment of inertia per unit length referred to wing 
elastic axis 

semi span of wing 

aerodynamic wing lift coefficient due to bending oscillations 
of the wing (see reference 5) 

aerodynamic wing lift coefficient due to torsional oscil- 
lations of the wing about its quarter chord (see reference 5) 


aerodynamic moment coefficient about wing quarter-chord point 
due to bending oscillations of wing (see reference 5) 

aerodynamic moment coefficient about the wing quarter-chord 
point due to torsional oscillations of wing about quarter- 
chord (see reference 5) 

mass per unit length 

static moment per unit length referred to wing elastic axis, 
positive for center of gravity behind elastic axis 

flutter speed, feet per second 

experimental flutter speed for wing without concentrated 
weight, 334 feet per second 

reduced flutter speed 

spanwise coordinate measured from wing root 
torsional component of ith coupled mode of vibration 
angular frequency at flutter, radians per second 
natural angular frequency of vibration in ith coupled mode 
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^1 + ig^j, where i is imaginary quantity \J - 1 


0h > 0h > 
1 2 

0 , 0 
V a 2 


uncoupled modal functions in first tending, second tending, 
first torsion, and second torsion, respectively 


mass density of air 


Subscripts: 

i, j designation of number of coupled modes; specific values 1, 

2, and 3 used for a particular coupled mode 


FLUTTER ANALYSIS WITH TEE USE OF COUPLED MODES 


The procedure for conducting a flutter analysis of the Rayleigh 
type for a given wing-weight configuration involves the selection of a 
set of modal functions to approximate the flutter mode, the formation 
of the flutter determinant, and the solution of this determinant for 
the flutter condition. 

The modal functions usually employed are either the coupled or 
uncoupled modes of vibration of the system. The term "uncoupled mode", 
as employed' in the present paper, refers to an imagined constrained mode 
in which, for pure bending, the chordwlse distribution of mass is con- 
sidered to act at the elastic axis of .the wing with no torsional defor- 
mation occurring. For pure torsion, the elastic axis Is considered 
restrained against bending. The term "coupled mode", as employed 
herein, refers to a combination of bending and torsional deflections 
appropriate to the natural harmonic vibrations of the freely oscillating 
undamped system. 

For the purpose of the present analysis, coupled modes have been 
selected. These coupled modes may be determined in any of a number of 
ways (see, for example, appendix II of reference 6). 

Once the coupled modes of the system are found, they are used 
together with the inertial characteristics of the system and the appro- 
priate aerodynamic coefficients to form the flutter determinant. The 
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complex flutter determinant, for the case of three coupled modes (see 
reference 6) , may he given in the following form: 



where the A^’s and C^j’s are generalized constants which are com- 
puted from the inertial properties of the system, the coupled modes, 
and the aerodynamic coefficients and are given by: 



The value is the angular frequency of the ith coupled mode- of 

vibration. The parameter H is a characteristic value given, in terms 
of the flutter frequency co and the concept of the structural damping 
coefficient g, by the relation 



where i, in this expression, is the imaginary quantity \J - 1. 

The functions h^ and refer to the bending and torsional 

components, respectively, of the ith coupled mode. 
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The flutter condition is determined from the nontrivial solution 
of this determinant. This solution may he obtained by various methods 
(see, for example, reference 7)- 


APPLICATION AND DISCUSSION OF RESULTS 


The method of analysis based on coupled modes discussed in the 
preceding section hq,s been applied to the case in reference 2 where the 
wing weight was designated as Ja. This case represented a uniform, 
unswept, cantilever wing, 48 inches long, with a concentrated weight 
mounted at various spanwise positions and is the configuration for 
which an analysis based on uncoupled modes (reference 4) gave the least 
satisfactory agreement with experiment. .The mass of the weight was of 
the same order as that of the wing. The position of the center of 
gravity of the weight was near the leading edge, well forward of the 
wing elastic axis. The calculations for flutter have been made for the 
wing without a weight and for the weight mounted at four different span- 
wise positions. 

The coupled modes of vibration employed in the present investi- 
gation were obtained by a process of matrix iteration based on computed 
influence coefficients. The wing with distributed mass was considered 
as a system of concentrated masses located at 12 equally spaced stations 
along its span. Where a concentrated mass was included in the system, 
its effects were considered at the appropriate spanwise station. The 
procedure used was essentially that outlined in appendix II of refer- 
ence 6. The actual iterative process was carried out on the Bell Tele- 
phone Laboratories X-66744 relay computer at the Langley Laboratory. 

In the solution of the flutter determinant the structural damping 
coefficients were assumed equal so that 

g-L = §2 = S3 = e 

where g is considered as a variable In the solution of the flutter 
determinant. The method of solution employed was that of reference 7* 

In this method the flutter determinant is put in the form of a set of 
simultaneous equations which are solved by a process of iteration. 

This iteration process yields values of g and cu from which the con- 
ventional plot of g against v can be obtained. For the cases con- 
sidered herein, flutter conditions for the case g = 0 were obtained. 

Analyses were made by use of the first three coupled modes of 
vibration of the system and by use of the three possible combinations 
of two coupled modes (that Is, first and second, first and third, and 
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second and third) . Designation of the modes by numbers is on a fre- 
quency basis, with the first mode being that which occurs at the lowest 
frequency, and so forth. A comparison of experimental and calculated 
coupled and uncoupled frequencies is given in table I. Note that the 
first, second, and third coupled modes correspond primarily to uncoupled 
first bending, first torsion, and second bending, respectively, except 
in the case of the no-weight condition. The results of the analyses 
are compared with the calculated results of reference 4 and the experi- 
mental results of reference 2 in table II and figure 1. Where "no solu- 
tion" is indicated in the table, either none exists dr it is well 
beyond the range of practical significance. In agreement with refer- 
ence 4 two uncoupled modes denotes uncoupled first bending and first 
torsion; three uncoupled modes denotes uncoupled first bending, first 
torsion, and second bending; and four uncoupled modes denotes uncoupled 
first bending, first torsion, second bending, and second torsion. 

For the wing without a weight, good agreement between calculated 
and experimental results was obtained with both coupled and uncoupled 
modes. For this case the computed results were not particularly 
dependent on the number of modes considered. As shown in reference 2, 
the elastic axis of the wing was very near the center of gravity, so 
that for the wing without a weight very little mass coupling existed. 

With the weight at the 11-inch spanwise station, the result 
obtained with two uncoupled modes (v = 1. 108 v Q ) iB slightly better 
than that obtained with the three coupled modes (v - 1.150 v Q ). Solu- 
tions were obtained for two of the three possible combinations of two 
coupled modes. When the first and second coupled modes (corresponding 
to uncoupled first bending and first torsion) were considered, a higher 
result (v = 1.450 v Q ) was obtained. A solution was also obtained when 
the first and third coupled modes (corresponding to uncoupled first and 
second bending) were considered. The result for this analysis was very 
high ^v = 2.075 v 0 ^. All of these results were well above experiment 
(v= 0.970 v Q ). 

With the weight at the 17- inch spanwise station, the nearest 
approach to the experimental result (v = 1. 144 v Q ) was obtained with 
the analysis based on four uncoupled modes (v = 1.491 v Q ). ® ie use 
of the three coupled modes gave a higher answer (v = 1.979 v 0 ) which 

was, however, nearer the experimental value than the result obtained 
with three uncoupled modes (v = 2.093 v Q ). No solution was obtained 
when only two uncoupled modes were considered. At this station a 
rather unexpected result was obtained in that an analysis based on two 
coupled modes (corresponding to the uncoupled first and second bending 
modes) gave a result (v = 1.680 v Q ) which was in better agreement with 
experiment than that obtained with three coupled modes. The signifi- 
cance of this result is not clear at present. 
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For the weight at the 46- inch spanwise station, the analytical 
result based on three coupled modes (v = 1.120 v Q ) and the result of 
the analysis based on the second and third coupled modes (v = 1. 078 v 0 ) 
were both in close agreement with experiment (v = 1.102 v Q ). For this 
case, the second and third coupled modes correspond to the uncoupled 
second bending and first torsion modes. Both of these results were 
considerably better than those of the uncouple d-mode analyses. When 
two uncoupled modes were considered, no solution was obtained. Con- 
sideration of three and four uncoupled modes both gave solutions, with 
that for three (v ='1.260 v ) being slightly higher than that for four 
(v = 1.228 v Q ). 0) 

With the weight at the 48-inch spanwise station (tip) the analysis 
based on three coupled modes gave a result (v = 1.072 v Q j in close agree- 
ment with the results of the analyses based on three and four uncoupled 
modes (v = 1.060 v Q and v = 1. 063 v Q , respectively), all of these being 
above the experimental result (v = 0. 958 v Q ). When the second and third 
coupled modes (corresponding to uncoupled second bending and first tor- 
sion) were considered, a result (v = 1.018 v Q ) in slightly better agree- 
ment with experiment was obtained. For this station, too, no solution 
was obtained from the analysis based on uncoupled first bending and 
first torsion. 

With very few exceptions analysis of these cases by either the 
coupled- or the uncoupled-mode approach gave results which were high in 
comparison with experiment. 

Figure 1 and table II indicate that increasing the number of modes 
in the coupled-mode analyses of these cases did not usually cause the 
result to converge toward the experimental results. In the uncoupled- 
mode approach of reference 4 the addition of a mode generally caused 
the result to converge toward experiment. 

When economy of computing time is considered, the uncoupled-mode 
approach is by far the better method. As an example, if only the deter- 
mination of the modes is considered, the computation of three coupled 
modes by the automatic computing methods employed required approximately 
23 hours and would correspond to a minimum of "JO hours of manual com- 
puting. In contrast to this, the time required to compute manually 
three uncoupled modes would be approximately 5 hours. The amount of 
time and labor required in the formation of the flutter determinant for 
the coupled-mode approach is also greater than that for the uncoupled- 
mode approach. 

For the cases studied, - the use of coupled modes yielded, in general, 
no better agreement with experiment than did the use of uncoupled modes 



2 


NACA TN 2375 


9 


and, in some cases, the uncoup led-mode approach was better. This result 
is in conflict with the frequently expressed thought that coupled modes, 
being more realistic insofar as the ground-vibration modes are concerned, 
should give a better approximation to the flutter mode and, therefore, 
a better result in flutter calculations. It should be recognized that 
the calculations performed were conducted for a uniform, unswept, canti- 
lever wing and that the conclusions may not necessarily be applicable 
to a sweptback wing. The observation that uncoupled modes give as good 
results as and, in some cases, better results than coupled modes is 
quite interesting. The reason for this result, however, is not known 
but it may be caused by the allowance of more freedom for the uncoupled- 
mode analysis to combine the bending and torsion modes. In the coupled- 
mode analysis, a mode has both a bending and a torsional component and 
these components are in phase and have fixed relative amplitudes. This 
phase and amplitude relation may cause some restriction in the combining 
of the modes and may be detrimental. 


CONCHJSIONS 


An investigation of the flutter characteristics of a uniform, 
unswept, cantilever wing of high aspect ratio and under conditions of 
high mass coupling has been made by means of an analysis of the Rayleigh 
type based on coupled modal functions. From the comparison presented 
herein of the results of the coupled-mode analysis with those of the 
uncouple d-mode analysis of MCA TN 1902 and with experiment, the fol- 
lowing conclusions can be drawn (which may not necessarily be applicable 
to a sweptback wing) : 

1. In most of the cases considered, which were selected to make a 
rather severe test of the use of modal functions, the use of either the 
coupled- or the uncoupled-mode approach gave results which were. high 

in comparison with experiment. 

2. Increasing the number of modes in the coupled-mode analyses of 
these cases did not usually cause the result to converge toward the 
experimental results. Such convergence with added uncoupled modes was 
indicated in MCA TN 1902. 

3. In comparison with the coupled-mode approach, the uncoupled- 
mode approach was by far the better method when economy of computing 
time is considered. 
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•4. For the cases treated herein, the use of coupled modes yielded, 
in general, no better agreement -with experiment than did the use of 
uncoupled modes and, in some cases, the uncoupled-mode approach vas 
better. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., March 22, 1951 
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TABLE I.- EXPERIMENTAL AND CAKHJIATED EREQUMCIES 


Weight 

• 

Span 

position 
(in. freon 
root) 

Frequency 

(epa) 

First mode 
(a) 

Second mode 
(b) 

Third mode 
(c) 

Experimental 

Uncoupled 

Coupled 

Experimental 

Uncoupled 

Coupled 

Experimental 

Uncoupled 

Coupled 

None 

0 

6.25 

6.55 

6.63 

35.3 

41.7 

41.2 

44.6 

48.5 

48.6 

7a 

11 

6.20 

6.53 

6.34 

22.8 

28.1 

25.4 

34.1 


39.0 

7a 

17 

5-90 

6.26 

6.21 

19.7 

23.0 

20.7 

(a) 

26.9 

35.5 

7a 

46 

3.16 

3.16 

3-17 

17.6 

l4. 7 

18.5 

31.8 

33.4 

33.6 

7a 

48 

2.70 

3-09- 

3.06 

18.0 

14. 5 

18.2 

29.8 

31.5 

31.9 


a ]Primary component of deflection similar to first cantilever bending mode. 

^Primary component of deflection similar to second cantilever bending mode in no~veight 
condition but similar to first cantilever torsional mode for all other veight conditions. 


c primary component of deflection similar to first cantilever torsional mode in no- veight 
condition but similar to second cantilever bending mode for all other veight conditions. 

^ot clear. 
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TABLE rr.- EXPiKDGHTAL AKD CALOTIATED EESDITS 


(a) Uncoupled modes. 



• 

Experimental 
results 
(reference 2) 

Calculated results for uncoupled modes 
(reference 

Weight 

Span 

position 
(in. frcm 
root) 

Two modes 
0hj* 

Three modes 

Four modes 

V V V \ 



f 

(cps) 

v/bCD 

v/v 0 

(a) 

f 

(cps) 

v/b© 

t A 0 

(a) 

f 

(cpb) 

v/bm 

(a) 

f 

(cps) 

v/txD 

V /t o 

(a) 

Hone 

0 . 

22.1 

7.22 

1.000 

25-2 

6.10 

H 

* 

6 

23-9 

6.8 

1.018 


\ 

H9 

7a 

Ta 

ii 

17 

17. * 

1*16.3 
|_ 26.8 

8.88 

*11.04 

7.02 

• 970 
|l.l44 

21. 3 
(c) 

8.28 

(c) 







mm 

(c) 

14.21 

23-52 

2.093 

32.8 

7.27 

1.491 

7a 

46 

21.8 

8.09 

1.102 

(<=) 

(c) 

(c) 

18.00 

7.64 

1.260 

18.2 

7.37 

1.228 

7a 

46 

21.4 

7-1* 

.958 

(c) 

(c) 

(c) 

25-50 

6.65 

1.060 

24.8 

6.06 

1.063 


(b) Coupled modes. 







Calculated results for coupled modes 

Wel^rt 

Span 

position 

Experimental 

results 

T»o nodes 


[in, from 

















root) 




Iflt Pnfl 

1st 

and 3rd 

2nd and 3rd 






f 

(cps) 

■e/hro 

v Ao 

(a) 

f 

(cps) 

T/brn 

v/v 0 

(a) 

t 

(cps) 

v/bco 

Vv 0 

(a) 

f 

(cps) 

v/boo 

v/v c 

(a) 

(cps) 

v/ba> 

v /v o 

(a) 

Hone 

0 

22.1 

7.22 

1.000 

(c) 

(c) 

(t>) 

24.2 

6.67 

1.009 

(c) 

(c) 

(c) 

24.2 

6.76 

1.024 

7a 

11 

17.4 

8.88 

.970 

14.5 

15.93 

1.450 

35.6 

9.28 

2.075 

(c) 

(0) 

(c) 

18.9 

9.68 

1.150 

7a 

17 

|*16.3 

|_26.8 

b 11.04 

7.02' 

|l.l44 

(=) 

(c) 

(c) 

30.0 

8.93 

1.680 

(c) 

(c) 

(c) 

30.6 

10.30 

1.979 

7a 

46 

21.8 

8.09 

1.102 

(<=) 

(c) 

(c) 

(c) 

(c) 

(c) 

24.3 

7.06 

1.078 

25.0 

7.1* 

1.120 

7a 

48 

21.4 

7.1* 

.958 

(c) 

(c) 

(c) 

(c) 

(c) 

(c) 

23.6 

6.88 

1.018 

2^.7 

6.92 

1.072 


\ - 334 fps. 

b 

In reference 2 with the weight at IT Inches from the root section the wing 
appeared to diverge. However, the oscillograph records for this case showed two 
possible flatter points, one corresponding to a frequency of 1 6. 3 cps (rather than 
the value of 16.0 cps recorded, in reference 2) and another corresponding to a 
frequency of 26.8 cps. Only the first of these values is noted in reference 2. 

°No solution. 
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Root Distance along span, in. Tip 


Figure 1.- Comparison of calculated and experimental flutter speeds for 
a particular wing-weight Bystem (weight 7a of reference 2) . 
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